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Abstract 



4d compact lattice QED in the quenched approximation. In the strong coupling 
phase we find exact zero-modes, quite similar as in non-Abelian lattice QCD. Sub- 
sequently we make an attempt to identify responsible topological excitations of the 
U(l) lattice gauge theory. 

Key words: Topology, Abelian gauge theory 
PACS: 11.15.Ha, 05.45.Pq, 12.38Ge 



In recent years the spectrum of the Dirac operator in QCD-like theories has 
been related to random matrix theory (RMT), see [1] for a review. RMT 
reproduces the effective, finite- volume partition function results of Ref. [2]. 
What enters into the RMT description of the low-energy, finite-volume scal- 
ing behavior are universal symmetry properties of the Dirac operator and 
the topological charge and the value of the chiral condensate of the gauge 
configuration under consideration [3,4]. Via the Atiyah-Singer index theorem 
the topological charge is mapped on the number of fermionic zero-modes of 
the Dirac operator. The symmetry properties of the Dirac operator fall into 
three classes, corresponding to the chiral orthogonal, unitary, and symplectic 
ensembles [4]. This classification is intimately connected to the chiral sym- 
metry properties of fermions. A good non-perturbative regularization of QCD 
should therefore retain those chiral symmetry properties. Based on the over- 
lap formalism [5] , such a lattice regularization emerged only recently with the 
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massless Neuberger overlap- Dirac operator [6] given by 



(1) 



where ^^Hyj{—m) is the usual Wilson-Dirac operator on the lattice and e the 
sign function. The mass m has to be chosen to be positive and well above 
the critical mass for Wilson fermions but below the mass where the doublers 
become light on the lattice. 

Previously the low- lying spectrum of the operator (1) was calculated for var- 
ious QCD ensembles [7] and agreement with RMT predictions was observed. 
In this paper we expand such calculations to compact, quenched lattice QED 
with the U(l) lattice gauge theory Wilson action 

S{Ui}^(5Y,{^-cose^) . (2) 
p 



Here j3 = Ui = f/^,^ = exp(?6'^^^), with 9^^^ e (-tt, +7r] and 9^ = 9.^^^^ = 
Ox,n + Ox+fi,y - Ox+i>,iJ. - Ox,y for 1/ 7^ /i. For /? < /^c 1.01, this theory is in the 
confinement phase, exhibiting a mass gap and monopole excitations [8]. This 
phase exhibits chiral symmetry breaking and the chiral condensate is, accord- 
ing to the Banks-Casher relation [9], determined by the small eigenvalues of 
the Dirac operator. Properties of the chiral phase have been studied numer- 
ically in Refs. [10-12]. In the strong-coupling limit /9 — 0, chiral symmetry 
breaking follows rigorously from infrared bounds [13] and it has also been 
calculated explicitly [14]. For (5 > /3c, the theory is in the chirally symmetric 
Coulomb phase with a massless photon [15]. There are many interesting ques- 
tions concerning the order of the transition between the two phases and the 
possibility of a non-trivial continuum limit for /? — > /3 ~ [16]. 

Using the staggered Dirac operator, RMT predictions for U(l) lattice gauge 
theory were previously investigated. In both, the strong coupling as well as 
the Coulomb phase, the nearest neighbor spacing distribution of the complete 
eigenvalue spectrum was found to be in agreement with the Wigner surmise 
of the unitary RMT ensemble [17], indicating quantum chaos. In the strong 
coupling phase, below the Thouless energy, the small eigenvalues were observed 
to contribute to the chiral condensate similarly as for the SU(2) and SU(3) 
gauge groups [18], in agreement with the chiral unitary RMT ensemble. 

For non-Abelian gauge theories the topological charge is determined by the 
instantons [19] of the gauge configuration. The topological structure of the 
Abelian U(l) gauge theory is different. Therefore, our first question is whether 
the overlap-Dirac operator will exhibit exact zero-modes [20]. To answer it, we 
have analyzed configurations on lattices at /3 = 0.9 in the confined phase 
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and at [3 — 1.1 in the Coulomb phase. With an overrelaxation/hcatbath algo- 
rithm we generated 500 configurations per lattice of linear size L = 4, 6 and 8. 
After thermalization, the configurations were separated by 100 sweeps, with 
each sweep consisting of 3 overrelaxation updates of each hnk, followed by 
one heatbath update. On each configuration the lowest 12 eigenvalues of the 
overlap-Dirac operator (1) were calculated using the the Ritz functional algo- 
rithm [21] with the optimal rational approximation of Ref. [22] to e {Hyj{m)) 
and m set to 2.3. 

In the confined phase exact zero-modes of the operator (1) were indeed found 
and are compiled in Table 1. There tiq'' denotes the number of configurations 
on which we found a zero-mode of degeneracy is the number of config- 

urations with no zero-mode. The total number of zero-modes in all produced 
configurations (of a given lattice size) is 



n, 



tot 



n, 



(3) 



The highest degeneracy observed was i/ = 3. No zero- modes were found in the 
Coulomb phase. 
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Table 1 

Exact zero-modes found in 500 configurations per lattice size L. The column Uq, 
= 0, . . . , 3, gives the number of configurations with u zero-modes, nQ°* is the total 
number of zero-modes and the last column presents the zero-mode susceptibility. 
The last row gives the Q values of the two-sided Kolmogorov test [25] and shows 
consistency of the RMT distributions (5) with our data. 

For each topological sector u, chiral RMT predicts the distribution of the 
lowest non-zero eigenvalue Amin in terms of the rescaled variable 

z^EVX^i^, (4) 



where V is the volume of the system and S the infinite volume value of the 
chiral condensate (ipip) determined up to an overall wave function renormal- 
ization. For the U(l) gauge group the unitary ensemble applies and the RMT 
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predictions for the u — 0,1,2 probability densities Pmhi(-^) [23,24] 



Pmin(^) 



Z 

2 



e for 1/ = 0, 
l2{z)e-^'/^ for 
.[72(^)'-/i(^)/3(^)]e-^V^ for u^2. 



(5) 



The chiral condensate is related to the expectation value of the smallest eigen- 
value. For degeneracy u we have 



where the result is supposed to be independent of u and 



(6) 



(^^^^) = jdzzp^±{z). 





(7) 



Using analytical and numerical integration the expectation values {z'^^'>) are 
easily calculated. One finds 



(^(°)) = 1.772453851 , {z^^^) = 3.107798700 
and (z^^)) = 4.344018125 . 



(8) 
(9) 



Using these numbers, the values E^''^ follow from Eq. (6) and Table 2 collects, 
for the various lattice sizes and z/'s, the results which we obtain from our 
data. Within our limited statistical accuracy, and allowing for some finite-size 
effects, the different E^'^'' values are consistent with one another. 
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Table 2 

Estimates of the chiral condensate T,^'^\ where the subscript u denotes the topolog- 
ical sector in which the estimate is obtained. The last column gives the weighted 
average over all topological sectors. 

Next, we compare the analytical distributions p^^j^ (5) with our computation. 
For z/ = 0, 1 and 2, Fig. 1 displays the exact RMT probability densities p^l^ 
of Eq. (5) and the corresponding histograms from our data. Using the defi- 
nition (4) of the 2;-variable and for E the E^^^^ values of Table 2, all lattice 
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Fig. 1. The exact RMT probability densities (5) for the lowest non-zero eigenvalues 
are compared with histograms from our data. 

sizes are combined for each u and the comparison is parameter free. Figure 1 
shows that the histograms follow the shift of the RMT probability densities 
and their general shape. The high peak of the v — 2 histogram is interpreted as 
a statistical fluctuation and consistent with the low statistics we have for this 
case. This claim is made quantitative by applying the two-sided Kolmogorov 
test [25] to the discrepancy between the RMT and our empirical distributions, 
which leads to the satisfactory Q values in the last row of Table 1. The topo- 
logical structure of the U(l) gauge theory is not as extensively studied as the 
non-Abelian case. Nevertheless, with torons [26], monopolc solutions [8] and 
Dirac sheets [27] a number of topological objects are known for the U(l) gauge 
theory and their possible relationship to the zero-modes of the overlap-Dirac 
operator is discussed in the remainder of this paper. 

Torons. For a fixed link (x, (x) we average the angle O^.f^ over the perpendicular 
3d space to obtain 6^ , i = 0, . . . , L—\. We define the toron charge for direction 
fi as the sum, using periodic boundary conditions, 

C^^E^M-i, (10) 

where 9^ ^_-^ is the shorter distance between 9^ and 9^_^ on the circle, yielding 
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an integer quantity. Table 3 gives the distributions of the total toron charges 



or 



(11) 



for all our configurations at /3 = 0.9 in the confined phase. No non-zero toron 
charges are found at /? = 1.1 in the Coulomb phase. 
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Table 3 

Distribution of the total toron charge (11). 

At /3 = 0.9 we checked for correlations between toron charges g^r and q^or with 
the degeneracy of our exact zero-modes and found none within the limitations 
of our statistical errors. Table 4 includes the average number of torons versus 
1/ (i.e. the average of the absolute values of the charges listed in Table 3). From 
L = 4 to 6 there is an increase with lattice size, but no statistically significant 
dependence on u. 

Monopoles and Dirac sheets. The U(l) plaquette angles 9,j.^^y arc decomposed 
into the "physical" electromagnetic flux through the plaquette 9x,ij,y and a 
number m-r^^u,^ of Dirac strings passing through the plaquette 



where O^^^u G (— 7r,-|-7r]. We shall call plaquettcs with rrix^^i, ^ Dirac pla- 
quettes. The monopole charges in elementary 3d cubes arc defined as the net 
number of Dirac strings entering or exiting these cubes. The worldlines of these 
monopoles on the dual lattice are closed, either within the lattice volume or 
by the periodic boundary conditions. The dual integer valued plaquettes 

''^x,tiv — 2^ixupa'>T^x,pa (13) 



form Dirac sheets bounded by the worldlines of monopole-antimonopole pairs. 
Due to the periodic boundary conditions they can also be closed without the 
presence of any monopoles or antimonopoles. In the remainder we shall refer 
to these as Dirac sheets. A Dirac sheet in the x-y plane of the Ad lattice must 
contain at least L^L^ Dirac plaquettes. Of course, monopole loops can form 
"holes" in these Dirac sheets. So, in practice we defined the number of Dirac 
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sheets in the x-y plane as the closest integer to Ni^2/{LzLt) where A^i,2 is the 
total number of Dirac plaquettes in the x-y planes. 
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Table 4 

Average numbers of torons, monopoles and Dirac sheets per configuration with fixed 
zero-mode degeneracy v. 

In the confined phase we find a large number of monopoles whereas the typical 
Dirac sheet numbers are and ±1 (only on the 8^ lattices we found at = a 
single configuration with two Dirac sheets). The average numbers of monopoles 
and Dirac sheets do not correlate with the topological charge u (i.e. our zero- 
mode degeneracy) as the numbers of Table 4 show. One might observe a slight 
increase of the monopole number with u which has no statistical significance. 
The average number of Dirac sheets shows no lattice size dependence within 
the limits of our statistics. In the Coulomb phase we find few monopoles and 
no Dirac sheets. 

All these topological configurations have in common with the zero-mode de- 
generacy that they arc turned on at the phase transition from order (Coulomb) 
to disorder (confined). However, in the confined phase we found no detailed 
correlation between any of the topological phenomena and the zero-mode de- 
generacy of the overlap-Dirac operator. This might be related to the fact that 
the zero-mode susceptibility, depicted in the last column of table 1, decreases 
and may approach zero for L — > oo. 

In contrast to these statistical findings there is an interesting observation that 
ordered Dirac sheet configurations give rise to zero- modes of the overlap-Dirac 
operator, with the number of zero-modes being equal to the number of Dirac 
sheets [5,20]. One might understand this from the fact that a Dirac sheet is a 
2d gauge configuration that contains unit topological charge in the 2d sense, 
kept constant in the two orthogonal directions. One might even verify this 
coincidence [5] by cooling or smoothing of equilibrium gauge fields [28] . 
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120 full fields 




1 20 monopole fields 




1 20 photon fields 




Fig. 2. Decomposition of the spectrum p(A) (normalized to the number of eigenval- 
ues) of the staggered Dirac operator into a monopole part Pmon(A) and a photon 
part Pphot (A) , averaged over 500 configurations on a 4^ lattice. 

To investigate these correlations further we have, following Refs. [29,30,10], 
factorized our gauge configurations into monopole and photon fields in the 
following way 



x' 



VfX ■ 



(14) 
(15) 



Here d' acts on x\ the quantities mx^^v and 6x,fiv are defined in Eq. (12) and 



G 



E.4sin2(2f ) 



(16) 



is the lattice Coulomb propagator. One can show that Q^^ii = + ^S*/?* 
is up to a gauge transformation identical with the original 9x,ij, defined by 
Ux,^, = exp(i6'^,^). 

This is a different way to probe the topological content of a gauge theory 
and we computed so far the total density of eigenvalues, with / p{X)d\ = V. 
We found that the zero-modes lie solely in the monopole part of the gauge 
field and are completely absent in the photonic field. Using periodic boundary 
conditions in space and anti-periodic boundary conditions in time [10], this 
was seen both for the overlap-Dirac operator and for the quasi-zero-modes of 
the staggered Dirac operator, see Fig. 2. At the moment we arc accumulating 
statistics to perform a decomposition into different topological sectors and to 
obtain an analogous analysis as we did with the original U(l) field above. It is 
of further interest to study space-time correlations between different topologi- 
cal objects. We pose the question of the existence of local correlations between 
the topological charge density and the monopole density. We plan to calculate 
spatial correlations with the wave functions related to fixed zero-mode number 
u and hope to understand the confining mechanism of compact QED. 
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